arXiv:astro-ph/0503573vl 26 Mar 2005 


Astronomy & Astrophysics manuscript no. astroph' scuba' grb3'1 February 2, 2008 

(DOI: will be inserted by hand later) 


SCUBA sub-millimeter observations of gamma-ray bursters 

Ml. GRB 030329: The brightest sub-millimeter afterglow to date 

I.A. Smith^ R.P.J. Tilanus^, N. Tanvir^, R.A.M.J. Wijers^, R Vreeswijk^, E. Rol^, and C. Kouveliotou^’^ 

* Department of Physics and Astronomy, Rice University, 6100 South Main, MS-108, Houston, TX 77005-1892 USA 
e-mail: iansmith@rice.edu 

^ Joint Astronomy Centre, 660 N. Aohoku Place, Hilo, HI 96720 USA 

^ Centre for Astrophysics Research, University of Hertfordshire, College Lane, Hatfield, Herts ALIO 9AB, UK 
'* Astronomical Institute ‘Anton Pannekoek’, University of Amsterdam and Center for High-Energy Astrophysics, 

Rruislaan 403, 1098 SJ Amsterdam, The Netherlands 
^ European Southern Observatory, Alonso de Cordova 3107, Casilla 19001, Santiago 19, Chile 

* Department of Physics and Astronomy, University of Leicester, Leicester LEI 7RH, UK 

’ NASA Marshall Space Flight Center, SD-50, NSSTC, 320 SparLman Drive, Huntsville, AL 35805 USA 

* Universities Space Research Association 

Received ; accepted 

Abstract. We discuss our ongoing program of Target of Opportunity (ToO) sub-millimeter observations of gamma-ray bursts 
(GRBs) using the Sub-millimetre Common-User Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope (JCMT). 
In this paper, we present all the ToO observations of GRB 030329. This was by far the brightest sub-millimeter afterglow seen 
to date. The flux density at 850 pm was approximately constant up to a break that took place ~ 7 days after the burst. This was 
consistent with being a jet break. The 850 pm results agree with those at longer wavelengths that show a brighter flux ~ 7 days 
after the burst, right at the time of the break. No short-lived large-scale brightenings were detected in the sub-millimeter light 
curve. However, the 850 pm light curve may have had a drop ^16 days after the burst. The peak of the afterglow emission was 
at ~ 90 GHz in the days before the break in the light curve. A simple modeling is consistent with the spectral indices remaining 
the same as the afterglow evolved, with the breaks in the spectrum moving to longer wavelengths at later times and the flux at 
the peak falling. No significant sub-millimeter emission was detected from the host galaxy. 
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1. Introduction 

The discovery of localized transients in the error boxes of 
gamma-ray burst (GRB) sources has led to intense multi¬ 
wavelength campaigns that have revolutionized our under¬ 
standing of these sources. For reviews see Van Paradijs et al. 
( I2000t and Meszaros i2002i . 

Sub-millimeter observations form a key element of the 
multiwavelength observations of the afterglow. They provide 
“clean” measures of the source intensity, unaffected by scin¬ 
tillation and extinction. We have therefore been performing 
Target of Opportunity (ToO) sub-millimeter observations of 
GRB counterparts using the Sub-millimetre Common-User 
Bolometer Array (SCUBA) on the James Clerk Maxwell 
Telescope (JCMT) on Mauna Kea, Hawaii. 

Including the results shown here, SCUBA has so far per¬ 
formed ToO observations of 23 bursts (Smith et al. ll999ll2()00l 
[200n [20051 Bloom et al.|T998l Galama et al.[T993 Kulkarni 
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et al.rr^ Frail et al. 120001120021120031 Berger et al. 120001 
12003 al Yost et al. I2002t . A complete summary of the previ¬ 
ous SCUBA results and the motivations for performing sub¬ 
millimeter observations of the afterglows and host galaxies is 
given in Smith et al. (|2005ll. In this paper, we present all the 
SCUBA ToO observations of GRB 030329. 


2. Multiwavelength observations of GRB 030329 

The long (> 100 sec) GRB 030329 was detected by HETE- 
2 and several spacecraft in the Interplanetary Network at 
20030329.484 UT (Vanderspek et al. l2004t . The intense ion¬ 
izing flux from the burst produced a sudden disturbance to the 
Earth’s ionosphere (Schnoor et al. IMBl . 

The redshift was determined to be 0.1685 (Greiner et al. 
I2003a> . corresponding to a luminosity distance of ~ 810 Mpc. 
This was one of the closest GRBs localized to date. Although 
its intrinsic luminosity was relatively low, the fact that it was 
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nearby made it one of the brightest bursts and afterglows ever 
recorded. 

GRB 030329 was an X-ray rich burst, and a distinct, 
bright, soft X-ray component may have been present during 
the burst. The X-ray afterglow seen by the Rossi X-Ray Timing 
Explorer and XMM-Newton was bright (Marshall & Swank 
[20U^ Marshall et al. EOO^ Tiengo et al. lIi^llXMll . However, 
the X-ray afterglow light curve was sparsely sampled. It might 
extrapolate back to match the soft X-ray tail seen during the 
burst. There appeared to be a break in the decay of the X-ray 
light curve ~ 0.5 days after the burst, and there may have been 
a flattening of the X-ray light curve ~ 40 days after the burst. 

Optical observations simultaneous with the burst did not 
detect the source, placing a limit of V ~ 5.5 on the reverse 
shock emission (Torii et al. unns- However, a bright optical 
afterglow was found, with R ~ 12.4 at 67 minutes after the 
burst (e.g. Price et al. 120031 Sato et al. 120031 Uemura et al. 
[200^ Urata et al. l2004b . 

The optical afterglow was polarized at the 0.3 - 2.5% level, 
and the polarization was variable on time scales down to hours 
(Greiner et al. I2003bt . The radio afterglow was not signifi¬ 
cantly polarized (Finkelstein et al. 120041 Taylor et al. l2005t . 

The optical afterglow decay can broadly be described using 
a broken power law with a break ~ 3 - 8 days after the burst 
(Lipkin et al. 120041 . However, extensive observations of the op¬ 
tical afterglow light curve found a complex evolution of bumps 
and wiggles superimposed on the overall fading (e.g. Greiner et 
al. l2003bl Smith et al. l2003l Matheson et al. 120031 Bloom et al. 
120041 Lipkin et al. l2004t . These could have been due to struc¬ 
tures in the forward and/or reverse shocks, or repeated energy 
injection from the central engine (e.g. Sari & Meszaros 120001 
Granot et al. l2003t . or inhomogeneities in the ambient medium 
that the fireball expanded into (e.g. Berger et al. I2000> . The in¬ 
terpretation of the various features in the light curve is therefore 
complicated. 

The radio and millimeter afterglows were the brightest 
recorded for any burst to date (Berger et al. I2003bl Sheth et al. 
120031 Kuno et al. 120041 Taylor et al. 120041 Kohno et al. l2005t . 
The radio and millimeter light curves were initially approxi¬ 
mately flat or slowly rising. After this, they broke to a power 
law Fy oc with a steep decay of d ~ 2. This is consistent 
with a jet break (e.g. Sari et al. ll999t . The light curves for the 
observations from 43 to 250 GHz all broke at ~ 8 days after the 
burst. In the single broken power law model, this would suggest 
that the jet break took place at ~ 8 days after the burst rather 
than at ~ 3 days which is allowed from the optical observa¬ 
tions. The peak occutTed at increasingly later times for longer 
wavelengths, for example, it was at ~ 30 days at 4.86 GHz. 
These longer wavelengths were below the synchrotron self- 
absorption frequency (Va) and/or below the typical synchrotron 
frequency for the lowest energy electron in the power law (v^). 
Interstellar scintillation was an increasingly important factor 
below 15 GHz for the early observations when the source was 
compact. At late times, the radio light curve showed an achro¬ 
matic flattening which is consistent with the blast wave becom¬ 
ing transrelativistic (Frail et al. l2005> . 

A more complex model for the multiwavelength emission 
used breaks in the early optical afterglow decay at ~ 0.25 and 


~ 0.5 days (e.g. Sato et al. IMHl Burenin et al. I2003t . The first 
break could be due to the cooling break frequency (vc) passing 
through the optical band. The second break would be similar 
to the one seen in X-rays; since it is achromatic, this suggests 
it is a jet break. However, an additional source of optical emis¬ 
sion would then be required ~ 1 day after the burst. A pos¬ 
sible model for the multiwavelength afterglow then uses two 
jets with different opening angles (Berger et al. I2003bl Sheth 
et al. 120031 Lipkin et al. 120041 Tiengo et al. I2004t . The ini¬ 
tial part of the optical and X-ray afterglow would come from 
a natTow ultra-relativistic jet. The late X-ray afterglow and op¬ 
tical brightening would come from a wide, mildly relativistic 
jet, with a jet break ~ 10 days after the burst. The plateau seen 
at longer wavelengths before the break could then be a con¬ 
volution of the falling flux density (as f from the natTow 
angle jet and a rising flux density (as from the wide an¬ 
gle jet. Although this model describes the overall trends in the 
light curves, the multiple bumps and wiggles in the optical light 
curve still remain to be explained. 

Beginning ~ 6 days after the burst it became evident 
that there was an underlying component from the Type Ic su¬ 
pernova SN 2003dh. The broad spectral features indicated a 
large expansion velocity similar to the Type Ic hypernovae SN 
1998bw and SN 1997ef (Stanek et al. l2ofel Hjorth et al. lMil 
Kawabata et al. 120031 Matheson et al. 120031 Mazzali et al. 
120031 Kosugi et al. I2004> . This was direct spectroscopic evi¬ 
dence that at least some classical GRBs originate from core¬ 
collapse supernovae. The SN 2003dh observations also indi¬ 
cated that the explosion was aspherical. 

The host is believed to be a low to moderate metallicity 
starburst dwarf galaxy with V ~ 22.7 and R > 22.5 (Fruchter 
et al. 120031 Matheson et al. l2003> . 

3. SCUBA observing details 

SCUBA is the sub-millimeter continuum instrument for the 
JCMT (for a review see Holland et al. I1999t . The observing, 
calibration, and reduction techniques used for GRB 030329 
were the same as has been used on other bursts and are de¬ 
scribed in detail in Smith et al. ifT^IMin 120051 . 

SCUBA uses two atTays of bolometers to simultaneously 
observe the same region of sky. The arrays are optimized for 
operations at 450 and 850 pm. Photometry observations are 
made using a single pixel of the arrays. The other bolometers 
in the arrays are used to perform a good sky noise subtraction 
(Archibald et al. mm . During an observation, the secondary 
is chopped between the source and sky. The term “integration 
time” always refers to the “on-i-off” time, including the amount 
of time spent off-source. An 18 sec integration thus amounts 
to a 9 sec on-source observation time. A typical measurement 
consists of 50 integrations of 18 seconds; we refer to this as 
a “run.” Each observation of a source in general consists of 
several such runs, with focus, pointing, and calibration obser¬ 
vations in between. 

Version 1.6 of SURF* (Jenness & Lightfoot l2000t was used 
to combine the integrations, remove anomalous spikes, flatfield 

' SURF is distributed by the Starlink Project. 
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Table 1. SCUBA 850 jum (350 GHz) afterglow observations of 
GRB 030329. 


Burst 

SCUBA observing times 

Time since 

Integration 

TgSO 

Afterglow 850 pm 


Start 

Stop 

burst (days) 

time (sec) 


flux density (mJy) 

GRB 030329 

20030403.383 

20030403.414 

4.914 

1800 

0.623 

37.8 + 4.5 


20030404.369 

20030404.401 

5.901 

1800 

0.467 

27.6 ± 2.9 


20030405.269 

20030405.302 

6.801 

1800 

0.356 

33.9 + 2.8“ 


20030409.222 

20030409.266 

10.760 

1836 

0.433 

10.9 + 5.1 


20030415.223 

20030415.254 

16.754 

1800 

0.300 

-0.2+ 1.9 


20030416.339 

20030416.370 

17.870 

1800 

0.192 

-1.9+ 2.7 


“ Not corrected for small pointing error 


the array, and subtract the signal from the sky bolometers. The 
zenith sky opacity was determined using “skydips” in which 
the sky brightness temperature was measured at a range of ele¬ 
vations. This was used to apply atmospheric extinction correc¬ 
tions to the observed target fluxes. At least one standard cali¬ 
bration target was observed each night to determine the abso¬ 
lute flux of the GRB. 

A typical integration time of 2 hours gives an rms ~ 1.5 mJy 
at 850 fim. However, the sensitivity depends significantly on 
the weather and the elevation of the source; since our ToO ob¬ 
servations are done on short notice, sometimes these factors are 
less than ideal. Based on observed variations of the gain factor 
and signal levels we estimate typical systematic uncertainties 
in the absolute flux calibrations of 10% at 850 /im. In general 
the rms errors of the observations presented here are larger than 
this uncertainty. 

The pointing of the JCMT is checked several times during 
the night to ensure that it is reliable. The pointing accuracy 
is usually a few arcsec. However, an error in the track model 
used between 2000 August 25 and 2003 April 25 resulted in 
pointing errors that were non-negligible (shifts larger than 4") 
for targets with elevations above 60° and over small ranges of 
azimuth. Since the 850 bolometric pixel has a diffraction 
limited resolution of 14", the target remains well within the 
bolometric pixel. However, there will be an error in the flux 
that is measured. For GRB 030329, the only problem was with 
a pointing calibration observation of 3C 273 on 2003 April 5. 
However, this did not significantly affect the result for that day. 

4. SCUBA observations of GRB 030329 

Unfortunately, SCUBA was not on the telescope at the time of 
the burst and it was several days before the first SCUBA obser¬ 
vation could be made. The weather was not good for most of 
the observations leading to larger rms uncertainties than nor¬ 
mal. Despite these difficulties, GRB 030329 was still easily de¬ 
tected at 850 jum in the initial observations. 

Table ^ summarizes all the 850 //m SCUBA ToO obser¬ 
vations for GRB 030329. These expand on and supercede the 
results in Hoge et al. (Unnii. The start times are the times when 
the GRB was first observed on each day, and the stop times 
are when the last observation of the GRB was completed. The 
time since the burst uses the mid-point between these start and 
stop times (although it should be noted that focus, pointing. 


and calibration observations may have also taken place in be¬ 
tween the start and stop times); the time is the elapsed time 
in the Earth frame, not the rest frame of the host galaxy. The 
“integration time” for each observation is the “on-i-off” time; 
only half of this is spent on-source. The zenith optical depth at 
850 /tm (t 85 o) is given for the time of the observation. 

The ~ 30 mJy source initially detected for GRB 030329 
was by far the brightest SCUBA afterglow seen to date. For all 
the previous ToO observations performed by SCUBA, the 850 
jum flux density from the afterglow and/or reverse shock was 
< 10 mJy. 

4.1. SCUBA 850 lum light curve 

The 850 /um light curve is shown in Figure^ This is similar to 
the observations at longer wavelengths. The flux density is ap¬ 
proximately constant initially, before breaking with a rapid de¬ 
cay. There is no evidence for an underlying dusty host galaxy. 

The weighted mean for the first three SCUBA observations 
was 30.7 +1.8 mJy at 850 fim. This gives the constant dashed 
curve shown in Figure[2for the first week. 

Unfortunately, we do not have high quality observations 
during the decay phase, due to the terrible weather at that time. 
However, the lack of detection 17 days after the burst allows 
us to use simple qualitative models to constrain the break time 
and/or the decay slope. 

The short dashed curve in Figure^shows a break at 8 days 
after the burst and a decay of d = 2. This would predict a flux 
density of 7.0 mJy at 16.8 days after the burst. This is therefore 
inconsistent with our non-detection at the 3.7cr level. 

To maintain 6 = 2 for the whole decay, it would be neces¬ 
sary for the break to occur earlier. The solid curve in Figure^ 
places the break at 7 days. This is just consistent with our 3cr 
limit at 16.8 days. 

To have the jet break occur at later times, it would be neces¬ 
sary to have a faster decay. For example, the long dashed curve 
in Figure n uses a break at 8 days and 6 = 3. Such a steep 
decay (on average) has not been seen in the other longer wave¬ 
length observations of this burst. However, as shown in §4.2, it 
is possible that a break moved through the sub-millimeter band 
between 7 and 17 days after the burst, and this might account 
for a faster decay. Alternatively, Kohno et al. (I2005> found that 
the flux density at 93 GHz was approximately constant from 
12 to 17 days after the burst before a rapid drop took place be- 
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Fig. 1. The 850 /tm light curve for GRB 030329. The same data 
and models are plotted in the linear-linear (top) and log-log 
(bottom) versions, except that the last two data points are plot¬ 
ted as 3cr upper limits in the latter. The initial dashed curve is a 
constant with Fy - 30.7 mJy. This breaks to Fy oc r^, with the 
break at 8 days and 6 = 2 (short dashed), the break at 8 days 
and 6 = 3 (long dashed), and the break at 7 days and 6 = 2 
(solid). 


tween 17 and 18 days after the burst. It is therefore possible 
that a similar rapid drop took place ^16 days after the burst at 
850 yum. 

The SCUBA flux density at 850 yum before the break in the 
light curve is consistent with being constant. Sheth et al. 
found a small (< 20%) drop in the flux density at 100 GHz 
starting April 2 followed by a rise to its mean level on April 
5. This is similar to our observation on April 4 being lower 
than on April 3 and 5. On April 5, the flux density at 43.3 GHz 
was the highest recorded, and the flux density at the four other 
longer wavelengths measured was also temporarily higher on 
that day (Berger et al. I2003bt . Thus this variability may have 
been real. If so, it complicates the determination of the actual 
time for the break. 

The observations during each SCUBA run show no signif¬ 
icant variability. One possible explanation for the fluctuations 
in the optical light curve is that the central engine had mul¬ 
tiple ejections, rather than a single episode (e.g. Granot et al. 
120031 . This would lead to refreshed shocks when the faster 
moving ejecta caught up with slower moving material. A pre¬ 
diction of this model is that there should be short-lived reverse 
shocks when the ejecta collide producing significant longer 
wavelength flares. The lack of large-scale variability at 850 yum 
agrees with the results from other observations (Sheth et al. 
I2003> . and argues against this model. However, we caution that 
our data only sparsely sample the light curve. 



Fig. 2. The longer wavelength spectrum for GRB 030329 be¬ 
fore the break in the light curve. The data points show typi¬ 
cal averaged values for ~ 6 days after the burst. Added to the 
850 //m (350 GHz) and 450 yum (670 GHz) data from this paper 
are the results at 4.86, 8.46 and 15 GHz (Berger et al. l2003bt . 
23 and 43 GHz (Berger et al. l2003bl Kuno et al. EITOi . 90 GHz 
(Kuno et al. 120^ . 100 and 250 GHz (Sheth et al. lMGll . The 
upper limit at 450 yum is 3cr. The solid curve is Fy = 0.15 
up to 15 GHz, Fy = 13.7 up to a peak of 61 mJy at 90 
GHz, and Fy = 582 v 02 peak. 

4.2. Evolution of the longer wavelength spectrum 

Due to the poor observing conditions, the source was not de¬ 
tected at 450 fim in any of the individual observations. The 
source was also not detected at 450 yum when the results from 
separate days are combined. The weighted mean for the first 
three days was -69 + 33 mJy at 450 /urn. 

Figure|2illustrates the longer wavelength spectrum before 
the break in the light curve using the average fluxes for the 
observations ~ 6 days after the burst. Since data has been 
combined from different observatories, and the source may be 
variable, the points will have uncertainties of up to ~ 20%. 
However, this does give a good qualitative picture for the longer 
wavelength spectrum at the time before the break. 

The solid curve in Figure |2] shows one qualitative fit to the 
data. The rise uses Fy oc up to 15 GHz and Fy cc y'U up to 
a peak of 61 mJy at 90 GHz. The fall is Fy oc v The short 
dashed curve extrapolates the curve to highlight the break 
at longer wavelengths. 

Spectral indices such as these are commonly found in syn¬ 
chrotron fireball models (e.g. Sari et al. Piran 

Wijers & Galama 119991 Chevalier & Li 120001 Granot et al. 
120001 Granot & Sari 120021 Panaitescu & Kumar E()04t . For ex¬ 
ample, the breaks could be = 15 GHz and v,„ = 90 GHz. 

For the R band (4.2 x 10^ GHz) this fit would extrapolate 
to give 0.9 mJy or R ~ 16. The source was that bright ~ 1 day 
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Fig. 3. The longer wavelength spectrum for GRB 030329 after 
the break in the light curve. The data points show typical av¬ 
eraged values for ~ 17 days after the burst. Added to the 3cr 
850 pm (350 GHz) and 450 pm (670 GHz) limits from this 
paper are the results at 4.86, 8.46, 15, and 43 GHz (Berger et 
al. l2003bt . 23 GHz (Kuno et al. EITOl . 93 GHz (Kohno et al. 
[2005}, 100 and 250 GHz (Sheth et al. l200^ . The solid curve is 
Fy = 0.294 up to 8 GHz, Fy = 9.41 up to a peak of 27 

mJy at 24 GHz, and Fy = 133 after the peak. The short 
dashed curve shows the fit from Figure |2 

after the burst, and was at F ~ 18 at 6 days after the burst. 
In this case, it would then be necessary to have another break 
in the spectrum between the sub-millimeter and optical bands. 
For example, this could be the cooling break v^. 

Figure 0 illustrates the longer wavelength spectrum after 
the break in the light curve using the average fluxes for the 
observations ~ 17 days after the burst. Again, the uncertainties 
in the points may be up to ~ 20%. The solid curve in Figure [^ 
shows one qualitative fit to the data, with Fy oc up to 8 GHz, 
Fy oc up to a peak of 27 mJy at 24 GHz, and a fall with 

Fy OC 

The short dashed curve in Figure |3 shows the fit from 
Figure |2 Qualitatively, the spectral indices have stayed the 
same, with the two longer wavelength breaks (va and v„, in the 
example above) moving to longer wavelengths and the flux at 
the peak falling. 

However, the model predicts a flux density of 7.1 mJy at 
850 pm at ~ 17 days after the burst. This is inconsistent with 
our lack of detection at the 3.7cr level if we use the limit from 
April 15, and at the 4.4cr level if we combine our last two 
SCUBA observations. 

One way to reduce the disagreement between the model 
and the SCUBA observation at ~ 17 days after the burst is if 
the third break (v^ in the example above) passed through the 
sub-millimeter band shortly before our last observations. This 


Table 2. Host galaxy fluxes for GRB 030329 determined by 
combining the SCUBA 850 pm and 450 pm observations taken 
on the last two days. 


Burst 

Redshift 

Host flux density (mJy) 
850 pm 450 pm 

GRB 030329 

0.1685 

-0.8 ±1.6 17 ±14 


could also help to explain why the decay of the 850 pm light 
curve was faster than 6-2. However, the spectral break would 
need to be sharp to fully explain the SCUBA limit. Also, in the 
simple fireball models, Vc should remain constant after the jet 
break. Thus this would require a more complex model, such as 
having the magnetic energy change with time (e.g. Yost et al. 

An alternative explanation for the lack of detection at 
850 pm at ~ 17 days after the burst is that there could have 
been a rapid drop in the light curve, as discussed in §4.1. The 
drop must have taken place earlier than 16.8 days after the burst 
at 850 pm. Figurej^shows that the 250 GHz point was still on 
the V F2 curve at 16.3 days after the burst. The drop took place 
between 17 and 18 days after the burst at 93 GHz. This in¬ 
dicates that if the drop is real, it did not occur at exactly the 
same time for all wavelengths. Instead, it swept from shorter to 
longer wavelengths. 

4.3. SCUBA host galaxy limits 

The results in Tablej^show the weighted means for the last two 
days of SCUBA observations. Although a small amount of af¬ 
terglow emission may be present, this contamination should be 
small. There is no indication for any sub-millimeter emission 
from the host galaxy. 

It is not surprising that GRB 030329 is not a bright qui¬ 
escent sub-millimeter source given that (1) it is at a low red- 
shift, (2) the host is probably a dwarf galaxy, and (3) the X-ray 
and optical observations show relatively little excess absorp¬ 
tion from dust above the Galactic value (Tiengo et al. 120031 
Matheson et al. EnOD- However, given the short observing 
times and less than ideal observing conditions, the limits on the 
host galaxy fluxes are less restrictive than for some of the other 
bursts whose host galaxies have been observed by SCUBA 
(Barnard et al. IM)1 Berger et al. I2003al Tanvir et al. IMmI 
Smith et al. [2005}. 

The SCUBA host galaxy compilation of Tanvir et al. ( I2004> 
focused on hosts with 850 pm rms flux densities <1.4 mJy. 
Since our rms is slightly higher than this, and the redshift of 
GRB 030329 is < 2, the lack of detection of a sub-millimeter 
host galaxy does not shed any new light on the question of 
whether GRBs are closely linked to the most luminous dusty 
star-forming galaxies. 

5. Summary 

GRB 030329 was by far the brightest sub-millimeter afterglow 
seen to date. Despite our limited observations, we were still 
able to constrain the break and/or decay light curve. Assuming 
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that the decay is not steeper than 5-2, the break must have 
taken place ~ 7 days after the burst at 850 /rm. This jet break 
could have been at a slightly later time if another break (per¬ 
haps the cooling break) passed through the sub-millimeter band 
between 7 and 17 days after the burst, or if there was a drop in 
the 850 fim flux at $ 16 days after the burst. 

No short-lived large-scale brightenings were detected in the 
sub-millimeter light curve. The flux density at 850 /rm was con¬ 
sistent with being constant up to the jet break. However, the 
850 fim results also agree with those at longer wavelengths that 
show a slightly brighter flux ~ 7 days after the burst, right at 
the time of the break. 

The peak of the afterglow emission was at ~ 90 GHz in the 
days before the break in the light curve. The spectral indices 
were qualitatively given by +2 and -1-1/3 up to the peak and 
-1/2 beyond it. A simple modeling is consistent with the spec¬ 
tral indices remaining the same as the afterglow evolved, with 
the spectral breaks moving to longer wavelengths at later times 
and the flux at the peak falling. 

GRB 030329 did not have any significant sub-millimeter 
contribution from the host galaxy. 

The observations of GRB 030329 show that for bright af¬ 
terglows it is important to monitor the flux evolution carefully 
to look for fluctuations in the light curve and to determine the 
breaks. Observations of new bursts are continuing to produce 
surprises, and there is much left to learn about GRB afterglows 
and host galaxies. To obtain a complete picture of their nature 
will require the careful study of many bursts to expand our sam¬ 
ple. Sub-millimeter observations with a ~ mJy sensitivity are a 
key component to the multi-wavelength coverage. To this end, 
our program of ToO observations using SCUBA is ongoing. 

Acknowledgements. The James Clerk Maxwell Telescope is operated 
by The Joint Astronomy Centre on behalf of the Particle Physics and 
Astronomy Research Council of the United Kingdom, the Netherlands 
Organisation for Scientific Research, and the National Research 
Council of Canada. 

We thank the JCMT Director Gary Davis for authorizing the ToO 
observations. We are indebted to all the observers whose time was dis¬ 
placed by these observations, and acknowledge the dedicated efforts 
of the JCMT telescope operators Jim Hoge and Jonathan Kemp for 
their valuable assistance with the observations. 

We are grateful to Scott Barthelmy and Paul Butterworth for 
maintaining the GRB Coordinates Network (GCN), and to the other 
ground-based observers for the rapid dissemination of their burst re¬ 
sults. 

We acknowledge the data analysis facilities provided by the 
Starlink Project which is run by CCLRC on behalf of PPARC. 

The work at Rice University was supported in part by 
AFOSR/NSF grant number NSF AST-0123487. 

References 

Archibald, E. N., Jenness, T, Holland, W. S., et al. 2002, MNRAS, 

336, 1 

Barnard, V. E., Blain, A. W., Tanvir, N. R., et al. 2003, MNRAS, 338, 

1 

Berger, E., Sari, R., & Erail, D. A. 2000, ApJ, 545, 56 

Berger, E., Cowie, L. L., Kulkarni, S. R., et al. 2003a, ApJ, 588, 99 

Berger, E., Kulkarni, S. R., Pooley, G., et al. 2003b, Nature, 426, 154 


Bloom, J. S., Erail, D. A., Kulkarni, S. R., et al., 1998, ApJ, 508, L21 
Bloom, J. S., van Dokkum, P. G., Bailyn, C. D., et al. 2004, ApJ, 127, 
252 

Burenin, R. A., Sunyaev, R. A., Pavlinksy, M. N., et al. 2003, 
Astronomy Letters, 29, 573 
Chevalier, R. A., & Li, Z.-Y. 2000, ApJ, 536, 195 
Finkelstein, A. M., Ipatov, A. V., Gnedin, Y. N., et al. 2004, Astronomy 
Letters, 30, 368 

Frail, D. A., Berger, E., Galama, T., et al. 2000, ApJ, 538, L129 
Frail, D. A., Bertoldi, F, Moriarty-Schieven, G. H., et al. 2002, ApJ, 
565, 829 

Frail, D. A., Yost, S. A., Berger, E., et al. 2003, ApJ, 590, 992 
Frail, D. A., Soderberg, A. M., Kulkarni, S. R., et al. 2005, ApJ, 619, 
994 

Fruchter, A., Levan, A., Hook, R., et al. 2003, GCN 2243 
Galama, T. J., Briggs, M. S., Wijers, R. A. M. J., et al. 1999, Nature, 
398, 394 

Granot, J., & Sari, R. 2002, ApJ, 568, 820 
Granot, J., Piran, T, & Sari, R. 2000, ApJ, 534, L163 
Granot, J., Nakar, E., & Piran, T. 2003, Nature, 426, 138 
Greiner, J., Peimbert, M., Estaban, C., et al. 2003a, GCN 2020 
Greiner, J., Klose, S., Reinsch, K., et al. 2003b, Nature, 426, 157 
Hjorth, J., Sollerman, J., Mpller, R, et al. 2003, Nature, 423, 847 
Hoge, J. C., Meijerink, R., Tilanus, R. P. J., & Smith, I. A. 2003, GCN 
2088 

Holland, W. S., Robson, E. L, Gear, W. K., et al. 1999, MNRAS, 303, 
659 

Jenness, T., & Lightfoot, J. F. 2000, Starlink User Note 216, Starlink 
Project, CLRC 

Kawabata, K. S., Deng, J., Wang, L., et al. 2003, ApJ, 593, L19 
Kohno, K., Tosaki, T, Okuda, T, et al. 2005, PASJ, 57, 147 
Kosugi, G., Mizumoto, Y, Kawai, N., et al. 2004, PASJ, 56, 61 
Kulkarni, S. R., Frail, D. A., Sari, R., et al. 1999, ApJ, 522, L97 
Kuno, N., Sato, N., Nakanishi, H., et al. 2004, PASJ, 56, LI 
Lipkin, Y. M., Ofek, E. O., Gal-Yam, A., et al. 2004, ApJ, 606, 381 
Marshall, E. E. & Swank, J. H. 2003, GCN 1996 
Marshall, E. E., Markwardt, C., & Swank, J. H. 2003, GCN 2052 
Matheson, T., Garnavich, P. M., Stanek, K. Z., et al. 2003, ApJ, 599, 
394 

Mazzali, P. A., Deng, J., Tominaga, N., et al. 2003, ApJ, 599, L95 
Meszaros, P. 2002, ARA&A, 40, 137 
Panaitescu, A., & Kumar, P. 2004, MNRAS, 350, 213 
Piran, T. 1999, Physics Reports, 314, 575 

Price, P. A., Fox, D. W., Kulkarni, S. R., et al. 2003, Nature, 423, 844 
Sari, R., & Meszaros, P. 2000, ApJ, 535, L33 
Sari, R., Piran, T., & Narayan, R. 1998, ApJ, 497, L17 
Sari, R., Piran, T, & Halpern, J. P. 1999, ApJ, 519, L17 
Sato, R., Kawai, N., Suzuki, M., et al. 2003, ApJ, 599, L9 
Schnoor, P. W., Welch, D. L., Fishm a n, G. J., Price, A. 2003, GCN 
2176 

Sheth, K., Frail, D. A., White, S., et al. 2003, ApJ, 595, L33 
Smith, D. A., Rykoff, E. S., Akerlof, C. W., et al. 2003, ApJ, 596, 
L151 

Smith, 1. A., Tilanus, R. P. J., Van Paradijs, J., et al. 1999, A&A, 347, 
92 

Smith, 1. A., Van Paradijs, J., Tilanus, R. P. J., et al. 2000, in Gamma- 
Ray Bursts: 5th Huntsville Symposium, ed. R. M. Kippen, R. S. 
Mallozzi, & G. J. Fishman (New York: AIP), 326 
Smith, 1. A., Tilanus, R. P. J., Wijers, R. A. M. J., et al. 2001, A&A, 
380, 81 

Smith, 1. A., Tilanus, R. P. J., Tanvir, N., et al. 2005, A&A, submitted 
Stanek, K. Z., Matheson, T., Garnavich, P. M., et al. 2003, ApJ, 591, 
L17 



I.A. Smith et al.: SCUBA sub-millimeter observations of gamma-ray bursters 


7 


Tanvir, N. R., Barnard, V. E., Blain, A. W., et al. 2004, MNRAS, 352, 
1073 

Taylor, G. B., Frail, D. A., Berger, E., & Kulkami, S. R. 2004, ApJ, 
609, LI 

Taylor, G. B., Momjian, E., Pihlstrom, Y., Ghosh, T., & Salter, C. 

2005, ApJ, in press lastro-ph/0412483 i 
Tiengo, A., Mereghetti, S., Ghisellini, G., et al. 2003, A&A, 409, 983 
Tiengo, A., Mereghetti, S., Ghisellini, G., et al. 2004, A&A, 423, 861 
Torii, K., Kato, T., Yamaoka, H., et al. 2003, ApJ, 597, LlOl 
Uemura, M., Kato, T, Ishioka, R., et al. 2003, Nature, 423, 843 
Urata, Y, Miyata, T., Nishiura, S., et al. 2004, ApJ, 601, L17 
Vanderspek, R., Sakamoto, T., Barraud, C., et al. 2004, ApJ, 617, 1251 
Van Paradijs, J., Kouveliotou, C., & Wijers, R. A. M. J. 2000, 
ARA&A, 38, 379 

Wijers, R. A. M. J., & Galama, T. J. 1999, ApJ, 523, 177 
Yost, S. A., Frail, D. A., Harrison, F. A., et al. 2002, ApJ, 577, 155 
Yost, S. A., Harrison, F. A., Sari, R., & Frail, D. A. 2003, ApJ, 597, 
459 


